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INTRODUCTION 

Coronal heating theories can be classified as either wave-heating or current-heating theories. 
According to the wave-heating model, convective motions below the solar surface generate 
MHD waves, which then propagate outward and dissipate their energy in the corona. Ob- 
servations of line widths in the solar chromosphere seem to rule out the possibility that the 
corona is heated by acoustic- or slow mode MHD waves (Athay and White 1978), leaving 
Alfven- and fast mode waves waves as possible candidates of energy transport. According to 
the current-heating model, the sub-surface convective motions cause random displacements of 
the photospheric magnetic footpoints, leading to twisting and braiding of the coronal magnetic 
field. The field-aligned electric currents associated with these twists are subject to resistive 
dissipation. The current-heating model applies only to ’’closed” magnetic structures such as 
coronal loops, whereas the wave model applies to both open structures (coronal holes) and 
closed structures (active regions). Another difference between the two models is the time scale 
t of the photospheric motions: in wave-heating models, r is smaller than or equal to the res- 
onance time 2L/va of the loop, whereas in current-heating models t » 2 L/va (L is the loop 
length, va is the Alfven speed). 

The details of the mechanisms responsible for wave- and current dissipation are presently 
not well understood. The problem is that dissipative processes such as plasma resistivity 
and viscosity are unimportant on the large spatial scales of observable coronal structures: 
dissipation can occur only if there are strong gradients in the magnetic- and/or velocity field, 
with length scales of 1 km or less in the corona. A crucial problem in any theory of coronal 
heating is, therefore, to explain how such small-scale structures are generated. In the context 
of the wave-heating theory, phase-mixing of Alfven waves, due to density inhomogeneities in 
the solar corona, provides a way to produce small-scale structures (Heyvaerts and Priest 1983, 
Sakurai and Granik 1984, Steinolfson 1985). 

In the current-heating theory, which is the subject of the present paper, the magnetic energy 
associated with the braided magnetic field must be similarly transferred to smaller scales. The 
process by which this "cascade” of magnetic energy occurs is not well understood. Parker 
(1972, 1979, 1983, 1986) suggests that the process is due to an intrinsic nonequilibrium of 
magnetic fields: the equations of magnetostatic equilibrium seem to allow solutions only for 
certain restrictive cases, in which the vertical component of the vorticity in the photosphere is 
essentially a constant of motion (Parker 1986). Since the motions on the Sun do in general not 
have this nice property, the magnetic field cannot simply adjust to the slow, random motions 
applied at the photosphere, but is forced to evolve on the Alfv6n time scale L/va • Parker 
assumes this dynamical relaxation leads to the formation of discontinuities (current sheets), 
where .magnetic reconnection will occur until the topological constraints are satisfied. Hence, 
according to the nonequilibrium model the formation of small-scale structures is due to a 
relaxation process that takes place on the Alfv4n time scale. 

Recently, I proposed a somewhat different picture of the cascade process (van Ballegooijen 1985, 
1986, hereafter papers I and II). According to this model there are no special restrictions on the 
velocity fields that may be applied at the photospheric boundary: the magnetic field in a coronal 
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loop can evolve slowly through a series of equilibrium states, without the necessity for magnetic 
reconnection. Hence, in our opinion the dynamical relaxation process proposed by Parker does 
not occur. The necessary condition for equilibrium found by Parker (1972J was shown to result 
from an incorrect ordering of terms in Parker’s perturbation scheme (cf. appendix of paper 
I). Therefore, a firm mathematical basis for Parker’s concept of nonequilibrium seems to be 
lacking. 

In the absence of nonequilibrium, the magnetic field evolves through a series of equilibrium 
states. Hence the formation of small-scale structures such as current sheets cannot be due 
to the relaxation process by which magnetic equilibria are reached; if current sheet formation 
occurs in closed coronal structures, it must be a result of the quasi-static evolution process. We 
expect, therefore, that the time scale for current sheet formation is related to the time scale t 
of the photospheric motions. To test this hypothesis, we need to understand the properties of 
braided magnetic fields as they evolve quasi-statically in response to random motions applied 
at the photospheric boundary. This is a rather difficult 3-dimensional problem, and a general 
procedure for computing statistical quantities such as the magnetic power spectrum are not 
known. However, one can gain some insight into the nature of the cascade process from a 
simplified 2-dimensional problem, in which only a single plane transverse to the mean magnetic 
field is considered (cf. paper II). In the following this model is briefly discussed. 

STATISTICAL MODEL 


Consider an initially uniform field Bo = Bq z, extending between two flat boundary plates 
located at z = 0 and z = L. We assume that the field is perturbed by a random, incompressible 
motion in the boundary plates, characterized by a correlation length t and a correlation time 
r (we assume l L). Then the velocity field in the interior of the volume is given by: 

v = \v x (x,y,z,t),v y {x,y,z,t),0\, (!) 

and the magnetic field is approximately given by: 

B = B 0 [b x {x,y,z,t),b y {x,y,z,t),l], (2) 

where b x and b y denote the transverse field components (b x ,b y C 1). Assuming ideal MHD, 
the induction equation can be written as: 


db x dv x 

dt dz flz ’ 
db v dv 

~dt = H = ° v ’ 


(за) 

( зб ) 


where djdz is the spatial derivative along fieldlines, and d/dt is the co-moving time derivative. 


The basic idea of the model is to consider an arbitrary plane z = zq in the interior of the volume 
(0 < zo < L), and to consider the velocity v(x, y, zo,t) and the velocity gradient a x (x,y,zo,t) 
at this plane as the independent statistical variables. Integration of equation (3a) then yields 
the transverse field b x : 

T 

6,(R,T) = f a x (r(t),t) dt, (4) 

o 

where r (t) = [x(t),y(t),zo] is the path in the zo-plane that ends at position R on time T. The 
correlation between the 6 z -values at two different points Ri and R 2 in the zo-plane is: 

c 6 (ar,t) =< 6 I (Ri,r)6 I (R 2 ,r) > 

T T ,-x 

= // < a I (r,(«'),(')a,(r 2 (l"),«") > dt' dt". 1 ' 

0 0 
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By making suitable assumptions about the statistical properties of a x , one can evaluate the 
right hand side of equation (5), which yields the magnetic correlation function C b (AR,T). 
We omit here the details of the derivation, which is given in paper II. The magnetic power 
spectrum P b {k,T) as function of transverse wavenumber k is obtained by taking the Fourier 
transform of C b with respect to AR. 
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The results of the calculation are displayed in the above figure, which shows the magn etic 
power spectrum (per unit log k) as function of dimensionless wavenumber kt. The different 
curves correspond to times T/t b = 0.4, 0.8, 1.2, 1.6, 2.0, and 2.4, where t b is the "braiding” 
time defined by: 

e 


and where £, r and u are the correlation length, correlation time and r.m.s. velocity of the 
photospheric motions, respectively. Note that magnetic energy, injected into the system at 
wavenumber k ~ £ -1 , is rapidly transferred to larger k : the maximum wavenumber in the 
spectrum increases exponentially with time, implying a rapid cascade of magnetic energy to- 
wards smaller length scales. This cascade takes place on the time scale t b , which is determined 
entirely by the statistical properties of the photospheric motions. Note, that t b is somewhat 
larger that the correlation time r, since the fluid displacements over one correlation time are 
generally smaller than the correlation length (ur < £). 

The exponential increases of the maximum wavenumber derives from the fact that, in an 
incompressible random flow, the separation between closely neighboring fluid particles Tift) 
and rj (t) increases exponentially with time; the e-folding time is of order t b . Hence, tne 
correlation function < Oi(ri(t),t)a x (r 2 (t),f) > appearing on the right hand side of equation 
(5) vanishes for time differences (T — t) larger than t b ln(£/ Ai2), when the separation of the 
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particles becomes larger than the correlation length. As a result the correlation function 
Cb(AR,T) has a logaxitmic dependence on AR, with a sharp peak at A R = 0. The Fourier 
transform of Cb therefore has significant power at high wavenumbers. 

The electric current density, j z = ( c/4ir)[dB y /dx — dB x /dy], involves the derivatives of the 
transverse field, and therefore the power spectrum of current-density fluctuations is given by 
Pj(k,T) = k I 2 Pb{k,T). It can be shown that the integral of Pj{k, T) over wavenumber, which 
is equal to < j 2 (T) >= j 2 ma , increases exponentially with time: 


jrms{T) ~ exp [2 v/27T7-1 • 

4 7T Lf L £& 


( 7 ) 


This should be contrasted with the "free” magnetic energy, which increases only quadratically 
with time: 


< Bl{T) > 


87T 


B 2 2 u 4 r 2 T 2 
8t r £2 L 2 ' 


( 8 ) 


Since the cascade of energy towards smaller scales proceeds exponentially in time, magnetic 
diffusion and reconnection will become important after a time t 1 which depends logaritmically 
on the magnetic Reynolds number: 


t\ « 


In R m 
4 6 v^' 


( 9 ) 


Here R m is defined as: 


Rm — 


P_ 

Vh' 


( 10 ) 


where 17 is the magnetic diffusivity based on the classical (Ohmic) resistivity. For the sun, 
Rm ~ 10 10 , so that t\ « 1.5t&. This implies that only a small number of braids can be 
introduced into the system before reconnection becomes important. 


I suggest that for t~>t 1 a statistically stationary state develops, in which there is a continuous 

transfer of magnetic energy from the scale £ where the energy is put in, to the scale IR “*/ 2 

where the energy is dissipated. The dissipation rate Eh in this stationary regime can.be 

estimated as the time derivative of expression (8), evaluated at the time t\ when reconnection 

processes first become important; this yields: 


b. Bo 2u 2 r lnR m 
E " ~ to 3i^ ' 

i.e., the heating rate depends logaritmically of the Reynolds number. Note that Eh is propor- 
tional to the product u 2 r, which is directly related to the effective diffusion constant of the 
photospheric motions: 

D = ^u 2 rV^r. (12) 

Observations of the spreading of active regions over time scales of months indicate that D 
is in the range of 150 - 425 km 2 /s (DeVore et al. 1985). With Bo — 100 G and L = 10 5 * * * * 
km, parameters typical for large active regions, we find Eh ~ 5 x 10~ 5 erg/cm 3 /s, which 
corresponds to an energy flux of 2.5 x 10 5 erg/cm 2 /s at each footpoint. This energy flux is 
a factor 40 smaller than the observed radiative- and conductive losses in active regions (cf. 
Withbroe and Noyes 1977). There are a number of possible reasons for this discrepancy: 

1) Our assumption that Eh equals the energy input rate at t = 1 1 probably under-estimates the 
heating rate, since the input rate may continue to increase for some time after reconnection 
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first becomes important. To see whether this is the case, it is necessary to include magnetic 
diffusion in the above analysis. 

2) There may be small-scale photospheric motions with D > 1000 km 2 /s which have so far 
escaped detection. 

3) It is possible that random photospheric motions are not the primary cause of coronal heating 
in active regions; periodic motion associated with Alfven waves may be a more important 
source of energy. 

In summary, we suggest that the quasi-static evolution of coronal magnetic structures is char- 
acterized by a cascade of magnetic energy to smaller length scales. This cascade process takes 
place on a time scale fj, determined entirely by the photospheric motions. The Ohmic heating 
rate Eh in the statistically stationary state was estimated using observational data on the dif- 
fusivity of photospheric motions; Eh turned out to be too small by a factor 40 when compared 
with observed coronal energy losses. However, given the fact that our theoretical estimate is 
based on a rather uncertain extrapolation to the diffusive regime, current heating cannot be 
ruled out as a viable mechanism of coronal heating. 
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